Extraction of electric field in heavily irradiated silicon pixel sensors 
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A new method for the extraction of the electric field in the bulk of heavily irradiated silicon pixel sensors is 
presented. It is based on the measurement of the Lorentz deflection and mobility of electrons as a function of 
depth. The measurements were made at the CERN H2 beam line, with the beam at a shallow angle with respect 
to the pixel sensor surface. The extracted electric field is used to simulate the charge collection and the Lorentz 
deflection in the pixel sensor. The simulated charge collection and the Lorentz deflection is in good agreement 
with the measurements both for non-irradiated and irradiated up to 10^^ Ueq/cm^sensors. 
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1. Introduction 

The properties of the silicon sensors designed 
for the CMS pixel detector 1 will change dur- 
ing the LHC operation. The innermost barrel 
layer of the CMS pixel detector is expected to 



be exposed to a fluence^ of 3 x 10 



14 



j/cm^ 



per year at full luminosity. The irradiation dose 
will be a few times larger in the case of the LHC 
luminosity upgrade. The silicon sensors behav- 
ior will be determined by the radiation damage, 
which changes the electric field in the silicon bulk 
and introduces charge trapping. This will lead 
to a reduction of the collected charge [2]. The 
pixel detector will operate in a 4 T magnetic 
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field and charge carriers will be deflected by the 
Lorentz force, which enhances charge sharing be- 
tween pixels and improves the spatial resolution. 
However, the bias voltage will be increased be- 
cause of irradiation and the Lorentz deflection 
will be reduced The spatial resolution de- 
pends on the charge collection, track position, sig- 
nal, noise and the Lorentz angle, and is degraded 
by irradiation. Here we present measurements 
of charge collection and Lorentz deflection as a 
function of depth in the silicon bulk for heavily 
irradiated pixel sensors. A new method for the 
extraction of the electric field in the silicon bulk 
is proposed and validated with a simple simula- 
tion. 
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2. Sensors and the measurement technique 

The sensors under study were designed for the 
CMS pixel detector and based on the "n-on-n" 
concept 0]. The bulk material is diffusively- 
oxygenated float zone (DOFZ) n-type silicon of 
(111) orientation and a resistivity 2-5 kJl cm. 
The pixels are formed by p-spray isolated n+- 
implants, while the p-n junction is formed by a 
large p+-implant on the back side. The thick- 
ness of the sensor is 285 ixm and the pixel size 
is 125 X 125 i^m^. The sensors were irradiated 
at the CERN SPS with 24 GeV protons at room 
temperature without applying bias voltage and 
then stored at -20°C. The tests were carried out 
at the CERN H2 beam line with 150-225 GeV 
pions. The beam entered the pixel plane at a 
shallow angle a =15° and the 3 T magnetic field 
was parallel to the beam (see Fig. The re- 
lationship between the z position of the created 
charge carriers and the corresponding arrival po- 
sition at the pixel plane along the x axis is given 
by z = xtana. Since the charge is always inte- 
grated in the pixel area, the smaller angle a is 
used the more precise location of charge carriers 
origin along the z axis is probed. The position 
of the beam exit point "O" was reconstructed in 
the pixel coordinates system using the beam tele- 
scope inj- The beam telescope consisted of four 
modules, each containing two silicon strip sensors 
measuring the x and y coordinates. The strips 
had a pitch of 25 [im, readout pitch of 50 (a.m 
and the spatial resolution of each plane was about 
l|j,m. The pixel sensors were bump-bonded to the 
PSI30/AC30 chip 0, which read out all signals 
from the 22 x 32 pixel matrix. The pixel sensor 
was cooled by Peltier elements down to -10°C. 
Both pixel and beam telescope signals were digi- 
tized using VME-based ADC modules controlled 
by a DAQ software written in Lab View and Lab- 
Windows/CVI (National Instruments). The trig- 
ger was provided by a PIN diode of size 3x6 mm^ 
placed between the beam telescope planes before 
the pixel detector. 

The electrons and holes produced by particles 
crossing the pixel sensor drift toward the elec- 
trodes. In absence of magnetic field the electrons 
are collected along the segment OA (see Fig. 




Figure 1. The defiection measurement technique. 



The holes move to the opposite direction and, to- 
gether with the electrons, induce the net current 
on the pixels situated along OA. In presence of 
a magnetic field charge carriers are defiected by 
the Lorentz force and the resulting current is in- 
duced on the pixels along the segment OD. This 
measurement technique was developed in 7 and 
used to measure the averaged Lorentz angle, O^, 
by fitting the deflection OD with a straight line. 
As we will see in section 01 the segment OD is 
curved, because the Lorentz angle depends on the 
electric field, which changes over the depth. How- 
ever, the experimental technique described in [7] 
can be applied for measuring the Lorentz angle as 
a function of depth in the sensor bulk. The bot- 
tom part of Fig. ^ shows the definition of (3{x) 
for an infinitely small section of the segment OD. 
Knowing the beam incident angle a=15° and the 
deflection angle /3{x) the Lorentz angle at a cer- 
tain depth 0l(z) is given by 

^ , , tandixiz)) ^/ X , 

tanOi z = ^ " = tan p(x) tana. (1) 

tan a 
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Figure 2. Asymmetry as a function of position in 
the xy plane for a non-irradiated sensor. 



The angle a is known, therefore, the Lorentz 
angle is calculated in each point in depth with 
the tangent to the segment OD, i.e. with the 
tan/3(a;). The geometrical position of the seg- 
ment OD in the xy coordinates plane can be de- 
termined from the signal asymmetry of two neigh- 
boring pixels. The asymmetry at the {x,y) posi- 
tion is defined as 



A{x,y) 



p/2 - Qx,y+p/2) 
{Qx.,y-p/2 + Qx,y+p/2) 



(2) 



where Qx,y-p/2 and Qx,y+p/2 is the charge col- 
lected in pixel, whose center is located at {x, y — 
p/2) and {x,y+p/2), respectively, andp =125 ^im 
is the pixel size. The asymmetry averaged over all 
events in each (x, y) bin is shown in Fig. |21 The 
asymmetry plot was divided into slices along the 
X axis. The charge spread is approximated with 
the Gaussian function, therefore the asymmetry 
in the z-th slice located at Xi (e.g. represented by 
the solid line in Fig. O was fitted with the stan- 
dard normal cumulative distribution function of 



y 



A{x = Xi,y) = cx 



{v-Vi)/s 



-'"/^dt ~ c, (3) 



where the parameter yi corresponds to the zero 
asymmetry position along the y coordinate for i- 
th slice (see Fig.EJ, c and s are the constant and 
spread parameters of the fit. The set of points 
{yi,Xi) determines the segment OD. The tan/3 is 
determined from the derivative dy/dx. The slight 



rotation (the line OA can be rotated with respect 
to the pixel row) of the sensor in the xy plane was 
subtracted using the data without magnetic field. 
Each point Xi corresponds to a certain depth via 
the relation Zi = Xi tan a and the deflection yi can 
be expressed as a function of depth. The mea- 
sured points were fitted with a 5-th order poly- 
nomial function (see Fig.^. 



3. Charge collection and the Lorentz de- 
flection 

Measurements without magnetic field were per- 
formed to determine the signal distribution along 
the segment OA see Fig. Assuming the av- 
eraged energy loss along the particle track to 
be uniform, the average signal in a pixel along 
OA is proportional to the charge collection effi- 
ciency originating at a certain depth in the sili- 
con bulk. The average charge collected by a sin- 
gle pixel as a function of the pixel position along 
OA is shown in Fig. 13 for a sensor irradiated at 
6.7x 10^'' Uoq/cm^. One can see, that even at low 
bias voltage (100-200V) some charge is collected 
from the p^ side and the charge has a minimum 
in the middle of the sensor thickness. Most of the 
charge, however, comes from the region close to 
the pixel implant. This behavior can be explained 
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Figure 3. Average collected charge as a function 
of the distance to the exit point for a sensor irra- 
diated at 6.7x 10^'' Uoq/cm^ for different bias volt- 
ages. 
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by a non-linear electric field and by trapping of 
charge carriers. If the sensor is operated in a mag- 
netic field, the charge carriers are deflected by the 
Lorentz angle Ql- The signal is induced on the 
pixels along the segment OD (see Fig. Since 
the electric field in irradiated sensor is not linear, 
the segment between the points "O" and "D" is 
curved, and in each point its tangent determines 
the deflection angle /3. The deflection of collected 
charge along the y axis as a function of the depth 
in the silicon bulk is shown in Fig ^ for different 
fluences and bias voltages. It must be noticed. 




z(nm) 

Figure 4. Lorentz deflection of the charge in a 3 
T magnetic field as a function of the depth. 



that at the edges of the silicon bulk (z < 17|xm 
or z > 268|J.m) the measured deflection has high 
systematic uncertainties (not shown in Fig.^J) due 
to the geometrical distortions of the electric field 
lines and incorrect reconstruction of the deflec- 
tion curve OD. The errors due to the particles 
multiple scattering were eliminated by selecting 
only straight tracks reconstructed by the beam 
telescope. 

4. Electric field measurements 

4.1. Electric field strength across the sili- 
con bulk 

In case of small Lorentz angles the low mag- 
netic field approximation can be used to describe 
the charge carriers motion in crossed electric and 
magnetic fields. This approximation physically 



corresponds to the situation where charge carri- 
ers travel only for a small arc of the circular orbit 
before scattering moves them from this orbit into 
another. The general expression of the current 
density in presence of crossed electric and mag- 
netic fields can be found in 8 . The Lorentz angle 
can be obtained from the direction of the current 
density as 

tSineL = rhfJ-B, (4) 

where is the Hall factor, B is the magnetic 
field and fi is the drift mobility. The charge car- 
riers mobility is determined by the lattice and 
impurity scattering. At the temperatures around 
300 K and for the impurity concentration up to 
10^^ cm~'^ the major contribution to the mobil- 
ity is from lattice scattering. For the impurity 
concentration up to 10^^ cm~'^ the impurity scat- 
tering contribution to the total mobility is below 
one percent and the Hall factor changes less than 
one percent in the range of impurity concentra- 
tion from 10^^ cm~'^ to 10^^ cm~^ Assuming 
that the irradiation induced defects act as im- 
purity atoms, the influence of irradiation on the 
mobility and Hall factor can be neglected for all 
irradiated sensors used for the tests. Therefore 
the measured Lorentz angle can be used to cal- 
culate the mobility using the Eq. 21 Most of the 
signal is due to the electrons contribution due to 
their shorter collection time and due to speciflc 
shape of the effective potential which is confirmed 
by the simulation (see Fig. Hand Fig.|HI). More- 
over, despite the fact that holes are also deflected 
during their drift, they will arrive to the back- 
side which has the same effective potential used to 
calculate the induced current on pixels using the 
Ramo-Shockley theorem JH]- Therefore the sig- 
nal induced by holes will be on the line OA even 
in presence of magnetic fleld and will not disturb 
the measured line OD which is only due to the 
collected electrons. Using the measured 8l(z) 
(tan Ol is a derivative of the deflection shown in 
Fig. Q and using Eq. ^ the electron mobility as 
a function of the depth is given by 
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where r/j = 1.15 is the Hall factor for the electrons 
and Bx = (3cos(15°)) T is the projection of the 
magnetic field along the x axis. The measured 
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Figure 5. Measured electrons mobility as a func- 
tion of depth for different fluences and bias volt- 
ages. The shaded regions correspond to the depth 
values where the mobility has large systematic 
uncertainties. 



electron mobility as a function of the sensor depth 
is shown in Fig|5l for different fluences and bias 
voltages. 

Using empirical parameterization of the field 
dependence of the electron mobility TT" and , 
one can derive the electric field as a function of 
the depth 

1 1/7= 



E{z) = 



1 



(6) 



where fie is the measured electron mobility and 
fiQe (low electric field electron mobility), Ece and 
7e are known empirical parameters from |11| 
and |12| . The parametrized mobility agrees 
within 5% with the measured mobilities |12, . The 
relative error of the electric filed ve calculated us- 
ing Eq.Elis related to relative error of the mobility 
r,, as 



1 







1 - 


(mo) 





(7) 



For the expected mobility range the error of the 
electric field extraction method is between 5 and 



15%. Fig. shows the electric field obtained ne- 
glecting the electric field lines distortion close to 
the pixel implants. The measurement is restricted 
to the depth range 17 |xm< z < 268 \xm for the 
reasons explained in section |21 The errors shown 
in Fig. are attributed to statistical fiuctuation 
of the collected charge and do not include the er- 
ror of the electric field extraction method, which 
is below 15%. For the non-irradiated sensor the 
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Figure 6. Extracted electric field as a function of 
depth for the non-irradiated and irradiated silicon 
sensors at different bias voltages. The shaded re- 
gions correspond to the depth values where the 
electric field has large systematic uncertainties. 



electric field is close to the classical linear field of 
an abrupt p-n junction. For the heavily irradiated 
sensors the electric field has a double peak with a 
distinct minimum close to the middle of the bulk. 
The origin of the double-peak electric field is qual- 
itatively described in |13| . A two-trap model pro- 
ducing a doubly-peaked electric field was imple- 
mented in a detailed detector simulation and the 
simulated charge collection was found to be in 
good agreement with the measurements |14ll5j . 
By integrating the electric field over the depth 
one can determine the potential drop across the 
silicon bulk. The potential drop agrees with the 
applied bias voltage within 15% for all sensors. 
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4.2. Cross-check of the measured electric 
field 

In order to check the measured electric field a 
simulation of the signal induced in the pixels was 
performed. The particle crosses the silicon sensor 
with an angle a = 15° (see Fig^ and the energy 
loss is assumed to be uniformly distributed. Nei- 
ther energy loss fluctuation nor charge diffusion 
was taken into account. In this simulation the 
electric field lines are assumed to be perpendic- 
ular to the silicon sensor planes and the electric 
field value as function of the depth was taken from 
the measurement shown in Fig. El The time de- 
pendent induced current is calculated using the 
Shockley-Ramo theorem ^U] 

i{t) - Qh{t)G[z{t)]-v!,[z{t)] + 

Qe{t)G[z{t)] ■ ve[zit)], ^ ' 

where and Qe are the holes and electrons 
charge values deposited by the particle energy 
loss, respectively, G is the weighting field, Vh and 
Ve the holes and electrons drift velocities, respec- 
tively. 

charge (el.) 




Figure 7. Simulated signal induced by holes. 



The drift velocity is calculated using the mea- 
sured electric field and the electrons and holes are 
trapped during the drifting time according to the 
exponential law 

Qhit) = Qohe-*/^N Qeit) = Qoee-*/"% (9) 

where the fluence dependent trapping probabili- 
ties Tfi^^ and Te~^ are calculated assuming a lin- 



ear dependence on irradiation fluence. The pro- 
portionality coefficients are 4.2x10""'^^ cm^/ns 
for electrons and 6.1x10^^^ cm^/ns for holes re- 
spectively 16 . For the non-irradiated sensors the 
trapping probability is set to zero, as the col- 
lection time is in the order of few nanoseconds. 
The induced signal was calculated separately for 
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Figure 8. Simulated signal induced by electrons. 



holes (see Fig. [Tj) and electrons (see Fig. ^ tak- 
ing into account the Lorentz force. The total in- 
duced signal is shown in Fig. The contribution 
from holes is significant only at the region close 
to the pixel implant while the total induced cur- 
rent is dominated by electrons. Figure [TUl shows 

charge (el.) 




Figure 9. Simulated total induced signal along 
the particle track projection with a magnetic field 
of (3cos(15°)) T. 
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the measured and simulated deflection as func- 
tion of depth. The charge deflection predicted 
by the simulation reproduces the measurements 
well. The simulation was performed also with- 
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Figure 10. Measured (markers) and simulated 
(solid lines) deflection as a function of depth for 
different fluences and bias voltages. 



out magnetic field to compare the charge collec- 
tion efficiency with the measurement. The simu- 
lation reproduces the measured values very well 
(see Fig. . The small discrepancies are due to 
charge diffusion, energy deposit fluctuation, elec- 
tric fleld distortion between the implants which 
were not implemented in the simulation. 

5. Summary 

A new method to extract the electric field in ir- 
radiated silicon pixel sensors is proposed and val- 
idated with a simulation. The method is based 
on a precise measurement of the Lorcntz deflec- 
tion as a function of depth in the silicon sensor 
bulk. The extracted electric field is used in a sen- 
sor simulation which reproduces very well both 
the charge collection and the Lorentz defiection. 
The method uses the electric field dependency 
on the mobility. Therefore, the precision of the 
method is limited by the precision of this depen- 
dence. However the empirical dependence agrees 
with the measurements very well |12j and the er- 
ror on the electric field is estimated to be less 
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Figure 11. Charge collected along the particle 
track for a sensor irradiated at 6.7 x 10^'' Uoq/cm^. 
The solid line is the simulated charge, dashed line 
is the measured one. 



than 15%. 

The Lorentz angle is used for mobility calculation 
and the error of the Lorentz angle measurements 
is dominated by the statistical fluctuations of the 
collected charge. The influence of the irradiation 
on the mobility and the Hall factor can be ne- 
glected for fluences up to lO^^Ueq/cm^ because of 
the lattice scattering dominates over the scatter- 
ing on the irradiation induced defects. 
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